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ABSTRACT 

Using a simple model of photodissociated atomic hydrogen on a galactic scale, it is 
possible to derive total hydrogen volume densities. These densities, obtained through 
a combination of atomic hydrogen, far-ultraviolet and metallicity data, provide an 
independent probe of the combined atomic and molecular hydrogen gas in galactic 
disks. 

We present a new, flexible and fully automated procedure using this simple model. 
This automated method will allow us to take full advantage of a host of available data 
on galaxies in order to calculate total hydrogen volume densities of giant molecular 
clouds surrounding sites of recent star formation. So far this was only possible on a 
galaxy-by-galaxy basis using by-eye analysis of candidate photodissociation regions. 

We test the automated method by adopting different models for the dust-to-gas 
ratio and comparing the resulting densities for M74, including a new metallicity map of 
M74 produced by integral field spectroscopy. We test the procedure against previously 
published M83 volume densities based on the same method and find no significant 
differences. The range of total hydrogen volume densities obtained for M74 is approx- 
imately 5-700 cm -3 . Different dust-to-gas ratio models do not result in measurably 
different densities. 

The cloud densities presented here add M74 to the list of galaxies analyzed using 
the assumption of photodissociated atomic hydrogen occurring near sites of recent star 
formation and further solidify the method. For the first time, full metallicity maps were 
included in the analysis as opposed to metallicity gradients. The results will need to be 
compared to other tracers of the interstellar medium and photodissociation regions, 
such as CO and CI1, in order to test our basic assumptions, specifically, our assumption 
that the HI we detect originates in photodissociation regions. 

Key words: galaxies: individual (M74) - galaxies: ISM - ISM: clouds - ISM: molecules 
- Ultraviolet: galaxies - ISM: atoms 



1 INTRODUCTION 

Molecular hydrogen is a much-investigated substance, as it is 
generally believed to be the key ingredient to star formation. 
Recent publications postulate a near-linear relation between 
(dense) molecular hydrogen gas densities and star formation 
rates (e.g. Krumholz et al.|2011||Schruba et aL 2011 1, at least 
on galactic scales, fueling this view. 

Yet molecular hydrogen itself is not easy to observe di- 
rectly, especially outside of our own Galaxy. The closest al- 
ternative, carbonmonoxide (CO), is a tracer of dense molec- 
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ular gas that is commonly used to infer molecular hydrogen 
colum n densities directly. Using 13 CO is preferred (Dickman 
19781, but 12 CO emission is stronger and therefore more 



readily detected. It is generally assumed that these CO in- 
tensities can be converted to molecular hydrogen column 
densities directly through a constant multiplication factor 
(Bloemen et al. 19861, commonly known as the X-factor. 



However, doubts about this conversion factor being truly 
constant across galaxies have been voiced since about the 



same time (Maloney & Black 19881. For example, Israel 
( |1997[ ) shows a clear metallicity dependence of the X-factor. 



The 'dark molecular gas' ( jWolfire et al.|2010| ), envelopes 
of molecular hydrogen without detectable CO emission, may 
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be missed and unaccounted for in the X-factor. A recent es- 
timate from observations of the Perseus molecular cloud in 
our own Galaxy by Lee et ah] ( |2012 \ is that approximately 
30 percent of the molecular gas is not accompanied by (vis- 
ible) CO. Ionized carbon [Cn] might be expected to trace 
this more diffuse molecular gas. Stacey et al. ( 1991 1 con- 



cluded that this carbon is formed in photodissociated gas 
between the giant molecular cloud and the ionized gas. They 
claimed that CO emission originates mostly from molecular 
gas exposed to elevated UV fields, not from cold disk clouds. 
Therefore CO may trace a combination of mass and excita- 
tion. For example, new results from the Herschel observa- 
tory in Mookerjea et al. (2011[) show how the CO likely gets 



photodissociated in the low-metallicity environment of M33. 
These observations of, among others, [Cn] are part of an on- 
going effort to explain its occurrence relative to the molec- 
ular (CO) and atomic (HI) gas. Against this background, 
it remains important and relevant to look into and develop 
different ways of estimating hydrogen abundances. 

An alternative approach to locate molecular clouds is 
to consider clusters of young, hot OB stars and the giant 
molecular clouds that surround them. These clouds are ei- 
ther the remnants of the clouds out of which these OB stars 
formed, or were newly formed when the gas was swept up 



by expanding giant HII regions. Boulesteix et al. ( 1974 \ de- 



scribe the (ring-like) HII regions of M33, which tend to be 
located on the edge of HI clouds with typical shell densities 
of the order of 1-10 cm -3 . Atomic hydrogen disks tend to 
show large HI holes with diameters ranging from 40 pc to 1 
kpc, where the holes tend to correlate well with OB associ- 
ations (M33: Deul & den Hartog 1990). Supershells remain 



actively investigated to date (see e.g. Dawson et al. 2011 



investigating supershells in our own Galaxy). 

In this context, the atomic hydrogen associated to the 
OB star clusters should be a direct product of photodissoci- 
ation. Certainly the OB stars surrounded by HI supershells 
conjure up an image of gas being blown into bubbles with 
molecular clouds forming in the shells, to be partially pho- 
todissociated by the UV radiation of the young stars. For a 



general review of photodissociation regions (PDRs), see Hol- 
|lenbach fc Tielens] p999 ) . In M83, atomic hydrogen appears 



to be photodissociated on a large scale ( Allen et al.||1986 l 
This property can be used to infer the total hydrogen gas 
density including the molecular hydrogen gas. For a recent, 
galaxy-wide application of this view of PDR-produced HI, 
see 



Higdon et al. (20111 



This procedure of locating potential PDR sites around 
OB star clusters, that we refer to as the 'PDR method', was 
presented previously for M33, M81 and M83 (|Allen et al 



1997| |Heiner et al.||2008a|b| |2011| ); see also |Smith et aT 



(20001 for M101. The method was investigated for selec- 
tion effects, resolution and distance effects. Until now, the 
analysis was carried out mostly by hand. This makes the pro- 
cedure slow to apply, and also potentially introduces biases. 
Therefore, we want to take the first step towards applying 
this method more routinely to a large number of galaxies and 
make the method widely available. To this end we present an 
automated 'pipeline' to extract total hydrogen volume den- 
sity measurements from a far-ultraviolet image, an atomic 
hydrogen image, and a dust model or dust map of a galaxy. 
These spot measurements present total hydrogen gas densi- 
ties at the location of potential PDRs across galactic disks. 



The densities are independent of CO measurements and in- 
dependent of any assumptions about the presence of molec- 
ular gas. They also do not depend on measurements (or the 
presence) of warm dust. We prefer this method because of 
the straightforward physics behind it, namely the process of 
photodissociation by ultraviolet radiation regulated by the 
local metallicity. 

For this first paper presenting the automated 'pipeline', 
we selected M74 as our first target, because it is practically 
face-on, reducing the influence of projection effects. Another 
reason to select M74 is the availability of detailed metallicity 
data, useful to increase the accuracy of the PDR method. 
This is the first time a metallicity map has been used in 
the method. M74 is more distant than M33, M81 or M83, 
but the tell-tale combination of far-UV sources surrounded 
by patches of HI emission is present. Paper II will feature 
a larger sample of nearby galaxies analyzed with our PDR 
method, enabling a more global evaluation of the method 
and its results. 

This paper is arranged as follows. First we briefly review 
the basic PDR theory, followed by the presentation of our 
automated measurement method. Then we present a newly 
obtained metallicity map of M74 that we use to derive the 
dust-to-gas ratio. We give the results of the PDR method 
for M74 using various metallicity measurements and com- 
pare the total hydrogen volume densities to those previously 
obtained in M83 ( |Heiner et al.|2008b| hereafter H08b). 



2 DERIVING TOTAL HYDROGEN DENSITIES 
FROM ATOMIC HYDROGEN 

This work is based on the assumption that a large part of 
atomic hydrogen in the disks of spiral galaxies is produced in 
photodissociation regions, as opposed to being 'primordial' 
and predating recent episodes of star formation. We briefly 
summarize this procedure, that comprises the PDR method. 

First, we need to locate so-called 'candidate PDRs', re- 
gions in galaxies dominated by clusters of young, hot OB 
stars. We expect these clusters to be surrounded by rem- 
nants of their parent giant molecular clouds, or alternatively, 
by recently formed GMCs in the expanding boundary of gi- 
ant extragalactic HII regions (GEHRs) with OB stars at the 
center. These regions, with a size of up to several hundred 
parsec, we collectively call candidate PDRs. Under the influ- 
ence of the far-ultraviolet radiation from the central cluster 
of stars, photodissociated atomic hydrogen appears on the 
surface of nearby GMCs. 

This PDR-produced HI is identified as locally enhanced 
HI column densities with a general HI background compo- 
nent. In general, several of these HI patches can be found 
per candidate PDR. Assuming a spherical geometry, we can 
then calculate the incident UV flux Go on these HI patches. 
Together with the local dust-to-gas ratio (DGR) we can then 
calculate the balance of photodissociation and therefore the 
total hydrogen volume density at the position of the HI col- 
umn. 

The final equation describing the balance of pho- 



todissociation is derived from Sternberg (1988) and Allen 



(20041, using a simple, two- level molecular hydrogen, one- 



dimensional approach. This allows us to calculate the total 
hydrogen volume density from the local atomic hydrogen 
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column density, the incident UV flux and the local dust- 
to-gas ratio. It includes recent improvements suggested by 



Heaton (20091. The result is 



n = 106- 



Gn 



(S/So) 



exp 



/ N HI (S/5 ) 
\ 7.8 x 10 20 



- 1 



(1) 



where the factor (<5/<5o)~ ' is just sligh tly different from th e 
factor (5/5o)~ O B used previously (e.g. Heiner et al. 20111. 
<5/<5o is the dust-to-gas ratio (scaled to the solar neighbor- 
hood value). Here, n is the total hydrogen volume density. 
Go is the incident far-UV flux scaled to the solar neighbor- 
hood value, and Nhi is the neutral hydrogen column density 
in cm -2 . In the context of a simple spherical geometry, the 
measured HI column density can be assumed to be equiva- 
lent to the HI column from central OB source in the direction 
of the GMCs surrounding it. 

The difference between the equation used in earlier ap- 
plications of this method and the current one is below the 
estimated uncertainty for normal values of S/Sq. 



3 AUTOMATED IMPLEMENTATION OF THE 
PDR METHOD 

The purpose behind developing the pipeline we present here 
is to provide the astrophysical community with a transpar- 
ent, yet configurable tool to apply the PDR method. 

The way the pipeline works can be roughly divided into 
three main steps. Firstly (Section 3.1l, candidate PDRs are 
identified. Secondly, (Section 3.2 1, measurements are taken 
from each candidate PDR (fluxes, geometry, column densi- 
ties and dust-to-gas ratios). Finally, total hydrogen volume 
densities are computed using the simple photodissociation 
model described in Section |2] 



3.1 Source selection 

Our method uses regions of strong FUV emission to se- 
lect candidate PDRs. We use the clumpfind algorithm of 



Williams et al. ( 1994 1 on FUV imagery to find and mask re- 



gions of enhanced emission. Clumpfind works by calculating 
contour levels on the selected image. The pixels contained 
within each contour level are labeled as belonging to the 
same region. When a new contour level is calculated, the al- 
ready defined regions lying within the new region are joined 
to the latter. Similarly, when two or more regions become 
adjacent, all of them are incorporated into the largest one. 
This process continues until each pixel belongs to one re- 
gion. The final product of clumpfind is a new file in which 
each pixel of the original image is labeled with the number 
of the region it belongs to. Hereafter, we will call this file a 
mask. 

Unless told otherwise, clumpfind will use all the pix- 
els in the image. In practice, a certain range of values must 
be excluded to avoid including the diffuse FUV background 
emission. We define this threshold as a function of the aver- 
age noise background level, <7n, which can be defined man- 
ually as an input parameter or interactively calculated. In 
the latter case, the user is asked to select a region of uniform 
FUV emission (i.e., neither a strong peak nor a region com- 
pletely free of FUV emission). <jn is then estimated through 



and iterative process: Firstly, the median and standard devi- 
ation, a, of all the pixels in the region is calculated; secondly, 
those pixels with values above or below 3cr are rejected (thus 
removing potential cosmic rays or "hot" pixels); this pro- 
cedure is repeated until the median value changes by less 
than 1%. Once on is known, the lowest contour level is de- 
fined as a number of times the background noise (n). While 
choosing an n too low (for example, n = 1) can ensure that 
no sources of emission are overlooked, spurious regions of 
background emission will likely be included. On the other 
hand, too strict a threshold value (n > 5) may miss some of 
the weakest regions, thus biasing the sample of HII regions. 
Another important parameter is the number of contours to 
fit. We have found that a number too large (for example 
50 contour levels) increases the computation time while not 
yielding any improvement (in some cases even producing 
worse results) . The optimal value will depend on the galaxy 
but our experience indicates that 20 contour levels (the de- 
fault value for the pipeline) produces good results without 
dramatically increase computation time. 

Our tests using different galaxies indicate that the re- 
sults of clumpfind are fairly robust against different sets of 
parameters, with the number of contours being the most 
critical value to obtain good results over a whole galaxy. De- 
pending on the chosen background level, far-UV clumps are 
identified at a larger galactic radius at a cost of potentially 
spurious detections. A balance must be struck by the user, 
who will need to assess whether any obvious sources of UV 
emission were missed and whether this can be remedied by 
changing the number of contours and the background level. 
We do not aim to extract a complete sample of candidate 
PDRs, but rather a representative sample in terms of far- 
UV flux and spread in galactocentric radius. The purpose of 
this method is to produce a statistically significant catalog 
of candidate PDRs in a way that minimizes the involvement 
of the user in the process. 

The initial mask produced by clumpfind has two issues 
that need to be dealt with. Firstly, some regions are broken 
into several parts that are still part of the same OB cluster. 
Secondly, even if regions are separate clusters they can be 
close enough to cause confusion. These issues can be solved 
by applying the following rejection criteria. First, we con- 
sider the typical size of a GMC and, if two regions are closer 
than a certain value (typically ~ 500 pc), the weakest one 
is rejected or joined to the strongest (the behavior can be 
defined as an input parameter). The "strength" of a region 
can be determined from its estimated total flux or the value 
of the peak pixel (again, configurable via an input param- 
eter). Similarly, when two regions are adjacent, a contrast 
test is carried out with the weakest region being rejected 
or joined to the strongest one. We find that this method al- 
lows us to resolve the great majority of the regions identified 
by clumpfind. This provides a selection of candidate PDRs 
similar to the more or less 'by-eye' method of selecting them 
used in our previous publications (like H08b). 

Finally, the user may opt to bypass the source detec- 
tion with clumpfind altogether and directly supply a list of 
coordinates. These coordinates could be, for example, a list 
of HII regions or a hand-picked list. This also means that 
the list of regions found by the clumpfind algorithm can be 
edited by hand if so desired. In this manner two use cases are 
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distinguished that can be mixed, namely a 'blind' selection 
of candidate PDRs and a pre-defined sample set of objects. 

Once the final FUV and HI regions have been selected, 
the pipeline produces catalogs in both VOTablsQ and DS9^ 
formats. FITS files containing each individual region in the 
far-UV and 21-cm are also created. The DS9 region file in 
particular facilitates easy inspection of the coordinates of the 
candidate PDRs, which is highly recommended considering 
the nature of the previously outlined selection process. 



3.2 Measurements and calculations 

Far-ultraviolet fluxes are measured using the task phot from 
the package digiphot in IRAF^] For each region, phot deter- 
mines the peak position of the candidate PDR and performs 
a series of concentric aperture photometries. Average fluxes 
per ring are calculated for each radius and the optimal aper- 
ture is determined where the average flux reaches its first 
local minimum after it has dropped to (by default, but con- 
figurable) below 50% the peak flux. Regions for which the 
optimal aperture cannot be determined using this method 
are rejected, which means that sources confused over larger 
areas and sources with a poor source contrast are excluded. 
The local background level is calculated using the average 
flux of the concentric aperture at the first local minimum. 

Based on the optimal aperture and local background 
emission, the pipeline calculates the net flux in the region as 
the cumulative flux at the given optimal aperture minus the 
estimated background emission for that area. The relative 
error ctfuv is defined as a fraction of the net flux, set as a 
parameter in the configuration file. 

Formally, we consider local peaks in the HI column 
density close to clusters of OB stars as PDR-produced HI 
patches. The column densities of these patches can be used 



directly in our calculations. We use SExtractor (Bertin & 
Arnouts|1996 1 to detect these patches automatically. While 
SExtractor was developed for star/galaxy recognition, we 
find it quite suitable for locating HI patches. In order to 
avoid confusion with larger scale structures, SExtractior is 
applied for each region on a subset of the HI map defined 
around the FUV peak. The SExtractor results are easily 
read and incorporated using Harry Ferguson's sextutils]^] 

The HI patches are corrected for diffuse emission using 
an overall background level which needs to be specified as 
an input parameter. The uncertainty am is set to half the 
HI background level. 

Using the recorded locations of each HI patch combined 
with the locations of the central UV sources, the source- 
patch separations pni are calculated on the plane of the 
sky. With pni we can calculate the incident flux Go and 
the source contrast at each HI patch, correcting for a (fixed) 
foreground extinction. The dust-to-gas ratio is determined 



1 |http: //www. ivoa.net/cgi-bin /twiki/bin/view/IVOA/ "| 
IvoaVOTable 

2 |http77Vhea- www. harvard. edu/RD/ds9/ 

3 IRAK is distributed by the National Optical Astronomy Ob- 



servatory, which is operated by the Association of Universities 
for Research in Astronomy (AURA) under cooperative agreement 

with the National Science Foundation. 

4 http: //www. stsci . edu/~f erguson/sof tware/pygoodsdist/ 
pygoods/sextu tils ,py| 



at the location of the HI patch. The pipeline provides three 
options for this: Using a simple galaxy-wide constant value, 
or a metallicity slope and a corresponding galactocentric 
radius for each patch, or an oxygen abundance map where 
a value can be read directly from the map. This last option 
in particular is expected to become more useful with the 
increasing availability of maps obtained from integral field 
spectroscopy. Finally, the total hydrogen volume density at 
each HI patch is calculated using Equation [T] We select the 
results to have a source contrast of at least unity. 



3.3 Configuration file 

One of the main goals we pursued while developing this au- 
tomated method was to make sure that the tool was highly 
customizable. Although we strived to provide reasonable de- 
fault values, most of the pipeline control parameters can 
be adjusted via a configuration file. The most important of 
these parameters are the number of levels used by clumpfmd, 
the background level of the atomic hydrogen map and the 
minimum acceptable signal-to-noise and several physical pa- 
rameters which control the FUV rejection mechanism. By 
providing sensible default parameters we aim to ensure that 
acceptable results are obtained with the user mainly provid- 
ing the input maps. The parameters mainly influence the 
localization of candidate PDRs and HI patches (through 
SExtractor) but not the actual physics that happens in the 
pipeline. 



3.4 Obtaining the software 

The pipeline is implemented using the Python programming 
languag^] The sofware will be released officially with Paper 
II (Sanchez-Gallego et al., in prep.), but in the meantime 
it can be obtained from the authors directly. The authors 
would be grateful to receive any feedback or suggestions that 
help to improve future versions of the tool. 



4 DATA 

4.1 Basic data and archival images of M74 

The local radiation field impinging on our candidate PDRs is 
estimated from far-UV photometry of the publicly available 



GALEX image of M74 ( Gil de Paz et al. 20071. We iden- 



tify atomic hydrogen features (that we expect to be PDR- 
produced) from the THINGS HI image flWalter et alTp008 |. 
The GALEX image has an angular resolution of approxi- 
mately 4 arcsec, whereas the THINGS image has a reso- 
lution of 6 arcsec (using the image produced with robust 
weighting). We list basic properties of M74 in Table [l] The 
adopted distance to M74 influences the deprojected sizes of 
our candidate PDRs. It also affects the value of R25 in units 
of kpc, but this has no influence on our results since it is just 
a scale of galactocentric radius. We chose this distance for 



consistency with Warren et al. (20101, who in turn adopted 



the distance measurement from the catalog of neighboring 



galaxies by Karachentsev et al. (20041. The foreground FUV 



http : //www . python . org 
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extinction is derived from the visual extinction E(B-V) us- 



ing the expression from Gil de Paz et al. ( 2007 1 and used to 



correct the measured GALEX fluxes. It is assumed that the 
candidate PDRs have an approximately spherical geometry, 
in which internal extinction between the central UV source 
and the HI patches is equal to the internal extinction to- 
wards the observer. In this scenario the internal extinction 
corrections cancel out and will be 0. If, due to projection ef- 
fects or because the spherical assumption is inaccurate, the 
extinction is non-zero, we will underestimate the incident 
UV flux and as a consequence underestimate the hydrogen 
density as well. We adopt a 'plausibility limit' of 500 pc, 
beyond which HI patches are assumed not to be associated 
with the central UV source. This limit compares favorably 
to a neutral gas disk scale height of 700 pc as adopted by 
Herbert-Fort et al. (20101 in their study of the outer disk of 



M74. 



4.2 Metallicity 

The balance of photodissociation is very sensitive to the lo- 
cal dust-to-gas ratio (sometimes abbreviated as DGR). We 
use the metallicity 12 + log(0/H) to derive the dust-to-gas 
ratio scaled to the solar neighborhood value 5/ So- After Issa 
et al. (19901, we assume a linear relation between the two 



(more recently, see Galametz et al.||2011 1 . After adopting a 



solar oxygen abundance of 8.69 ( Allende Prieto et al.|2001[ ) a 
dust-to-gas ratio can be derived. Typically, metallicity gra- 
dients are available in the literature for spiral galaxies. Since 
the scatter in such measurements is considerable, and likely 
intrinsic ([Rosolowsky & Simon 20081, local metallicity val- 
ues are preferred. In the case of M74, we present a new 
metallicity map that we used to improve the total hydrogen 
density measurements. We also use the metallicity map of 
M74 published by Sanchez et al. (20111 and their derived 



metallicity gradient to compare the results. 



4-2.1 Metallicity slope 

Assuming a constant slope in the metallicity of M74 as a 
function of galactocentric radius, we can obtain a value of 
the oxygen abundance (and by extension the dust-to-gas 
ratio) for every candidate PDR. The most recent result for 
M74 comes from Rosales-Ortega et al. ( 2011 1, based on their 



M74 metallicity map. Of the various methods they describe, 
we adopt the one using the 03N2 index from Pet tini fc Pagel| 
(20041: 



\og(0/H)/dex P 25 = -0.66, 



(2) 



where P25 is the optical radius of M74 (in arcmin). 

After correcting for the different distance to M74 and 
using a solar oxygen abundance of 8.69 we arrive at 



log(J/*o) 



-0.060X.R + 0.24, 



(3) 



with R in kpc. By adopting the 03N2 result we maintain 
consistency with the metallicity map we present here that is 
derived using the strong line relations presented in the same 
paper. 



4-2.2 A new metallicity map of M74 

The M74 metallicity map presented in this paper was 
produced using a data cube obtained in September 2008 
with the imaging Fourier transform spectrometer SpIOMM 
I Grandmont et ar]|20 03 Dris sen et al.|[2008| |Bernier et al 



2008 Charlebois et al.|2010 1 attached to the 1.6 m telescope 



of the Observatoire du Mont-Megantic. A red filter, covering 
the wavelength range from 6480 to 6820 A and a folding or- 
der of 16 with 206 steps was used for a spectral resolution of 
around 5 A and a total exposure time of 3 hours 15 minutes. 

The SpIOMM data cube is composed of interferograms. 
An image of the galaxy is seen on each of them, and as a 
whole, they allow us to get a spectrum for each pixel of the 
image. The detector used was a Princeton Instrument EEV 
CCD camera with 1340 x 1300 pixels with a 12' circular field 
of view. The data were binned 2x2 during the readout in 
order to match the size of a pixel (1.1" when binned) with 
the seeing (of the order of 1.5"). This gave us a total of 
435500 spectra. 

The first step of the data reduction was to align the in- 
terferograms using the centroids of field stars. Then, for each 
interferogram, we subtracted a bias and a dark image (col- 
lected with the same instrument). We used a dome flat image 
(with the instrument mobile mirror fixed away from the zero 
path difference position) to calibrate the pixel-by-pixel re- 
sponse of the detector. Cosmic rays were subtracted in each 
interferogram. The sky background variation was corrected 
by considering an image region away from the galaxy and de- 
void of other objects. In this region, we measured the average 
flux value between two consecutive interferograms (one with 
a constructive and one with a destructive interference). All 
consecutive interferograms should have presented the same 
average sky value, but because of changing conditions dur- 
ing the long exposure, this was not the case. Therefore the 
interferograms were corrected in intensity to take into ac- 
count the measured variation seen in the sky. It should be 
noted that as of 2011, SpIOMM has been equipped with a 
second CCD camera, drastically improving sensitivity and 
permitting measuring the total flux values. This allows cor- 
rections for sky variation to be made during the observa- 
tions. Finally, the Fourier transform was applied to the data 
cube and the wavelength calibration was done using a high- 
resolution data cube of a reference HeNe laser at 6328 A. All 
these operations were performed with IDL (the Interactive 
Data Language). 

Line measurements for the HII regions should be done 
on the average spectrum for each regions in order to in- 
crease the signal-to-noise ratio. Therefore, HII regions have 
been identified using the Hllphot program of |Thilker et al.| 
( |2000 1, created specifically for this purpose. Using the SpI- 
OMM data, we first created an Ha map, a red continuum 
map, and a combined Ha+continuum map of the galaxy, 
where Hllphot calculated HII region contours. The rele- 
vant parameters we used for Hllphot were: Nllpercent = 
1% (percentage of [Nil] contamination in the Ha flux) and 
size_max = 500 pc (maximum size of the HII regions). A 
total of 541 HII regions have been identified. They range in 
size from 4 to 93 pixels. An example interferogram is shown 
in Figure [2] At this point, we did not directly combine the 
spectra included in each region, but reapplied the Fourier 
transform to the data cube after combining the pixels corre- 
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sponding to the HII regions in each interferogram. This step 
also improved the signal-to-noise ratio in the final spectrum 
for each HII region. The sky background level and lines are 
removed by building an average sky spectrum using all pix- 
els away from the galaxy and other objects. An example HII 
region spectrum is shown in Figure [3] 

Since the goal here was to obtain a metallicity map from 
measurements of the two nearby lines of Ha and [NIIJA6583, 
no flux calibration and extinction corrections have been done 
yet (a complete analysis of these data, along with a blue 
cube, will be part of another paper by Rousseau-Nepton). 
The two lines have been measured using a Gaussian fit with 
the IDL task gaussf it. The uncertainties for these measure- 
ments are dominated by the position of the continuum. 

The metallicity (Z) was finally calculated using the 
[Nil] /Ha ratio and the polynomial relations of |Pettini fc| 
|Pagel| ( |2004 1. Altough this ratio gives a good estimate of 
the metallicity, it is sensitive to the ionization factor (q). 
The degeneracy effect is even more important at high Z. We 
calculated the error due to this degeneracy effect using the 
relations between the [Nil] /Ha ratio and the Z, q value that 
can be found in |Kewley fc D opita (2002). Figure [I] shows 
our metallicity map. 

We compare our metallicity map to the one obtained 
with the PPAK IFS Nearby Galaxies Survey (PINGS), as 



this radius (Zaritsky & Christlein 20071. Additionally, the 



presented in Sanchez et al. (2011 1 and Rosales-Ortega et al 



(20111. There are a number of differences. While both the 



SpIOMM and PPAK instruments are integral field spectro- 
graphs, they are constructed quite differently. SpIOMM is a 
fourier transform imaging spectrograph with a field of view 
large enough to capture the entire optical extent of M74, 
but at a cost of wavelength coverage. PPAK is a fiber spec- 
trograph with a comparatively limited field of view and a 
mosaic of M74 was constructed with it, covering a smaller 
area than SpIOMM was able to cover. PPAK's higher sensi- 
tivity (mostly due to it being mounted on a bigger telescope) 
and wavelength coverage allowed various ways of analyzing 
the results. In this paper, we simply use the resulting oxygen 
abundance 12 + log(0/H) map, with 



Iog(<5/5 ) = (12 + \og{0/H)) - 8.69. 



(4) 



A detailed comparison between the two different metal- 
licity maps is beyond the scope of this paper. The full SpI- 
OMM M74 metallicity map and comparisons with other 
published data will be presented in a different publication 
(Rousseau-Nepton et al., in prep.). 



5 HYDROGEN DENSITIES IN M74 
CANDIDATE PDRS 

Here we present total hydrogen volume densities calculated 
as in Eq. [I] near candidate PDRs across the disk of M74 
using different sources for the local dust-to-gas ratio. 

The candidate PDRs are located in the manner de- 
scribed in Section [3] The selected regions are shown in Fig- 
ure [4] (the circles indicate the HI patches surrounding the 
candidate PDRs). We also created far-UV/HI overlay plots. 
An example is shown in Figure [5] We note that the regions 
found by clumpfind are all inside the optical radius R25 
due to our choice of a detection threshold, although M74 is 
known to have a young stellar population out to two times 



metallicity maps are limited to the optical radius as well, 
which is the focus of this work. 

Fig. [6] shows from top to bottom the measurements 
obtained from our pipeline: the far-UV fluxes of candidate 
PDRs across M74, the associated HI column densities found 
near the UV sources, the separation between the central 
UV sources and the associated HI patches and the resulting 
incident UV fluxes Go at each HI patch. The results were 
selected to have a source contrast Go/G{, 9 of at least 1. 

The full dataset is available electronically, but an exam- 
ple is shown in Table[2] The datasets are in VOtable format, 
and the candidate PDRs are numbered consistently in the 
different M74 files of the densities derived from the metal- 
licity slope, PINGS metallicities and SpIOMM metallicities. 
Not all candidate PDRs located through clumpfind result 
in reported total hydrogen volume densities, since they are 
dropped if the source contrast is too low, no HI patches 
are found nearby or if no metallicity data is available at the 
source location. For the PINGS map we did not have individ- 
ual uncertainties so they are not included. The uncertainties 
in the metallicity slope values are estimated conservatively 
at 10%. 

The measured HI column densities increase steadily 
with galactocentric radius before flattening out, although 
the exact trend is unclear. The HI column densities peak 
near 3 x 10 21 cm -2 , reflecting perhaps a beam smoothing 
issue similar to the one found in HI column densities ob- 
tained from M83 in H08b. A constant HI background level 
of 2.9 x 10 20 cm -2 was subtracted, which has very little 
influence on the final total hydrogen volume density. 

The distribution of the values of pni is not as discrete 
as in our previous work (e.g. Fig. 4 in Heiner et al.|2011 1 due 
to a different measurement technique. Instead of manually 
locating HI patches within increasing annuli at discrete sep- 
arations, we have SExtractor locate the HI patches and use 
their exact location to calculate the separation, deprojected 
using the galaxy's position angle and declination as before. 
The accuracy of the measurement of this separation is still 
limited by the angular resolution of the HI map and, more 
importantly, by the unknown 3-dimensional position of the 
HI patches with respect to the central UV source. 

The values of Go do not show a strong radial trend, 
which is consistent with previous observations of other 
galaxies. 

We compared three different sources of the dust-to-gas 
ratio (see Fig.[7|: the SpIOMM metallicity map (black dots), 
and the metallicity map and gradient from Rosales-Ortega 



et al. ( 2011[ ). The gradient is based on their metallicity map, 



so it follows the actual range of metallicities. It can be seen 
that the metallicities from the SpIOMM metallicity map are 
systematically lower across most of the range of galactocen- 
tric radii. 

The total hydrogen volume densities computed from the 
three different metallicity sources are shown in Fig. [8] No sig- 
nificant difference is present, although the values obtained 
from the SpIOMM metallicity map are consistently higher 
because of the lower dust-to-gas ratios. The median differ- 
ence in dust-to-gas ratios for candidate PDRs present in 
both metallicity maps is 0.38, where the PINGS dust-to-gas 
ratios are mostly higher. The median difference in total hy- 
drogen volume densities is about 10 cm -3 . The overall range 
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of densities is between ~ 5cm and ~ 700cm with a few 
higher density outliers. 

Finally the plot with the observational selection effects 
(Go vs. n) is shown in Fig. [9] (only for the results based on 
the SpIOMM metallicity map). This figure is similar to the 
Allen et al.| ( |2004[ ), e.g. their Fi gure 3, but the axes 



an issue. The full results on which the plots are based are 
available through the electronic version of this journal. 



Fig. 10 (left) shows the HI patches surrounding the can- 
didate PDRs located by clumpfind. We reproduce Fig. 1 of 
H08b showing the candidate PDRs identified in that paper 



are switched. As can be seen in their figures, constant HI 
column densities correspond to straight lines in this type of 
plot. The Roman numerals indicate the selection effects: 

I. This is the HI column density upper limit of 5 x 
10 21 cm -2 at a typical 5/5o of 0.6. Theoretically, the HI 
column becomes optically thick at 5 x 10 21 cm -2 , but this 
value is not observed due to beam smoothing (see H08b). If a 
value of 3 x 10 21 cm -2 were used, the line would be closer to 



the point cloud. As Allen et al. ( 2004 I noted, for typical ob- 



servations the optically thick limit is a few times 10 cm" 
depending on spin temperatures, HI profile linewidths and 
local opacities. We also show the position of this limit for 
three additional S/80 of 0.3, 0.5 and 0.7 (dash dotted lines), 
where the highest value corresponds to the lower-right line. 

II. The HI lower limit of 2.9 x 10 20 cm" 2 at a typical 5/S 
of 0.8 is a sensitivity limit. The dash dotted lines show the 
same limit using <5/<5o of 0.3, 0.5 and 0.7, where the lowest 
value corresponds to the upper-left line. 

III. The lowest volume density that we can observe as a 
consequence of the radio beam diameter and the HI sensi- 
tivity limit. For M74, it is approximately 0.4 cm -3 . This 
limit is independent of our method and in practice we don't 
obtain values this low due to the dust-to-gas ratio being low 
(relative to e.g. M83). 

IV. The minimum Go (~ 0.01) that we can use is a combi- 
nation of the lowest measured flux and the largest acceptable 
value of phi- In practice a stronger cut-off is the source con- 
trast. The values of Go drop below one at a value of about 
log(Go) = -1.2, explaining the actual boundary of the points. 
This is an observational bias, since the source contrast could 
be improved with more sensitive data. 

V. The maximum Go (~ 1400) that we would obtain from 
the data if the highest flux coincided with the smallest mea- 
sured separation. As with the previous limit, it is dependent 
on the observational data. A higher value could be achieved 
with higher resolution data, to measure smaller separations, 
or if higher UV fluxes were measured. Source contrast is not 
an issue in this situation. There are not many flux measure- 
ments towards the high end because sources this bright are 
rare, but additionally it is feasible that the highest fluxes 
cause HI columns too low to be detected. Therefore, higher 
resolution would need to go together with higher sensitivity 
in order to be able to obtain points towards this limit. 



6 BENCHMARKING AGAINST PREVIOUSLY 
OBTAINED M83 DENSITIES 

To better understand the performance of our automated 
density analysis, we compute hydrogen densities for M83 
using the same input data as were used in H08b (see Table 
[T|. Due to certain limitations in our code, we did not use a 
variable position angle like we did in that paper, but fixed 
it at 43 degrees. The results are also limited to the extent 
of the optical disk, making the fixed position angle less of 



as Fig. 10 (right). Note that in the former figure there are 
multiple circles per candidate PDR, corresponding to indi- 
vidual HI patches that are assumed to be PDR-produced, 
whereas in the latter figure the central UV sources, or can- 
didate PDRs, are marked. 

Results from the automated PDR method for M83 are 
shown in Fig. [TT] From top to bottom: Dust-to-gas ratio, 
measured far-UV flux, HI column densities, separations be- 
tween the HI patches and the central UV sources and the in- 
cident flux Go • The dust-to-gas ratio is determined directly 
from the metallicity gradient. The far-UV fluxes show the 
usual distribution, diminishing in strength going outwards. 
The HI column densities as a function of galactocentric ra- 
dius appear to peak around ~ 0.35 R25, which is slightly 
sooner than what was observed in H08b, although this can- 
not be determined very accurately due to the scatter in the 
column densities. The HI background column density was 
set to 2.0 x 10 20 cm -2 and subtracted from all measured 
HI patches, consistent with H08b. The separations between 
the central source and the HI patches cover the full range 
of distances (up to 500 pc). This limit was referred to as a 
'confusion limit' in H08b, where it was set to 22 arcsec, or 
almost 500 pc. The results of the automated method were 
also selected to have a source contrast of at least unity. This 
is slightly different from H08b, where no such selection was 
performed, but results with a source contrast below 0.5 were 
marked. In principle, the source contrast is a measure of how 
dominant the UV source is with respect to a general back- 
ground UV radiation field. A value of unity then means that 
the source is just as bright as the background. 

The total hydrogen volume densities obtained here are 
compared to the ones from H08b in Fig. |13| with the den- 
sities from this work in black, and the H08b results in gray 
(open and closed circles). Only the range of galactocentric 
radii for which measurements were obtained using the au- 
tomated PDR method is shown. No significant difference in 
the range of values can be seen, although on average fewer 
HI patches were identified by the automated method than 
by eye in H08b. Since we chose to use only densities from 
HI patches experiencing a source contrast above 1, this plot 
shows the H08b results differentiated by this threshold: the 
open circles are the H08b densities with a source contrast 
below one. There is no clear separation between the two dif- 
ferent sets of densities, although most of the lower-contrast 
densities have lower density values. Comparing the (loga- 
rithmic) densities of the two sets with a source contrast 
above 1, we find a median value of 1.25 and 1.15 for the 
H08b densities and the densities obtained here respectively. 
This corresponds to median densities of 17.6 and 14.1 cm -3 
respectively. Neither set of (logarithmic) densities exhibits 
a significant correlation coefficient (0.35 and -0.04 respec- 
tively), meaning that neither set can be shown to have a 
correlation between galactocentric radius and (logarithmic) 
total hydrogen volume densities. 

For completeness, we show how the observational se- 
lection effects limit the results in Fig. |12| with the roman 
numerals denoting the same limiting factors as in H08b. The 
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limits to the HI column densities are slightly different be- 
cause of the different range of galactocentric radii covered, 
and the typical 8/Sq values are different as a result as well. 

I. The HI column density upper limit of 3 x 10 21 cm -2 at 
a typical <5/(5o of 1.9 is appropriate for the range of values 
obtained here. Theoretically, the HI column becomes opti- 
cally thick at 5 x 10 21 cm -2 , but this value is not observed 
due to the resolution of the HI map. 

II. The HI lower limit of 3.8 x 10 20 cm~ 2 at a typical 
S/6o of 2.2 is a sensitivity limit. It is slightly higher than 
the limit adopted in H08b, since lower HI column densities 
were identified by eye than automatically by SExtractor. 
The typical 8 /So for this value is also different from the one 
in H08b because in the latter case it was taken at a much 
larger galactocentric radius. 

III. The lowest volume density that we can observe as a 
consequence of the radio beam diameter and the HI sensi- 
tivity limit. For M83, it is approximately 0.7 cm -3 . Due to a 
slightly higher adopted HI lower limit, this value is slightly 
higher than the one in H08b. 

IV. The minimum Go (~ 0.05) that we can use. It is 
higher than the value in H08b (0.005) due to the lowest mea- 
sured UV flux by the automated method being higher than 
in H08b. Those lower fluxes originate mostly from larger 
galactocentric radii (not used here). 

V. The maximum Go (~ 770) that we can obtain from 
the data. Obtaining this value would require a chance su- 
perposition of both the highest UV flux being accompanied 
by a very close HI patch. It is somewhat higher than the 
H08b value (565) due to a higher maximum measured flux 
(since we do not aim to obtain a complete sample of can- 
didate PDRs this can happen), but the actual value would 
only be obtained by chance. 



7 DISCUSSION AND CONCLUSIONS 

The main factor influencing the uncertainty in the result- 
ing total hydrogen volume density is the dust-to-gas ratio. 
In principle the use of local metallicity values instead of an 
averaged slope should improve the accuracy of the results, 
but this is countered by the intrinsic scatter of the metallici- 
ties of individual HII regions. In other words, the individual 
points should have a smaller error, but the overall result do 
not change. 

It should also be noted that the distance to M74 is far 
from certain. For example, the distance as quoted in |Rosales^| 
Ortega et al.l 12011b puts M74 a full 2 Mpc further away. 



However, the effects of this uncertainty on our results are 
limited: because we use a 'plausibility cut-off' when consid- 
ering HI patches surrounding candidate PDRs, adopting a 
larger distance would effectively lower this cut-off, reducing 
the number of HI patches under consideration. Other mea- 
surements are not affected. For example, the value of Go 
is dependent on the ratio between pni and the distance to 
M74. Finally, a different distance could influence the selec- 
tion of candidate PDRs (since we reject the ones that are 
too close together), potentially leading to a higher density 
of selected candidate PDRs. 

We reiterate that our method provides spot measure- 
ments of total hydrogen volume densities. Due to the nature 
of the one-dimensional, plane-parallel PDR model, we do not 



obtain information about cloud morphology, nor does the 
treatment provide a fraction of molecular hydrogen. Addi- 
tional assumptions about cloud size would have to be made 
that are beyond the scope of this work. What we obtain 
are total hydrogen volume densities of GMCs surrounding 
clusters of OB stars, probing the interstellar medium of the 
galactic disk wherever these clusters are present. When we 
compare the range of our input parameters and hydrogen 
volume densities to the theoretical perspective provided by 
Krumholz et al. (20081, it seems that our values are reason- 



able for clouds as they are found in our Galaxy. Additionally, 
we find no evidence of changing cloud densities as a function 
of galactocentric radius. 

One of the most fundamental assumptions in our ap- 
proach is taking the HI patches near clusters of OB stars as 
produced by photodissociation and, more specifically, that 
these observed HI columns are connected to the nearby 
young stars. With our new automated method we can gather 
the kind of multi-galaxy statistics that will allow us to test 
this underlying assumption by comparisons with tracers of 
star formation and photodissociation regions. 

The findings presented in this paper can be summarized 
as follows: 

• We present an automated procedure to calculate total 
hydrogen volume densities in nearby spiral galaxies, based 
on the 'PDR method'. This procedure can replace the pre- 
vious manual application of the PDR method. 

• We apply this method using metallicity maps for the 
first time. As compared to using metallicity slopes, the den- 
sities resulting from using metallicity maps are more accu- 
rate but not significantly different from the ones obtained 
using a slope. This implies that in general the densities ob- 
tained from the PDR method can be expected to be similar 
regardless of whether a metallicity slope or map is used, 
which is key to applying this method to larger samples of 
galaxies. 

• We show a new metallicity map of M74 obtained with 
the SpIOMM instrument at the Mont-Megantic telescope 
and compare it to a similar map from the literature. The 
SpIOMM map is limited to the [Nil] /Ha ratio, but fully 
covers the extent of the M74 optical disk. 

• The range of total hydrogen volume densities obtained 
in M74 is within 5 - 700 cm -3 . 

• We apply the method to M83 and compare the results 
to the earlier results of H08b in order to validate the auto- 
mated method. Within the more limited range of galacto- 
centric radii, the results are not significantly different. We 
are therefore confident that we can proceed to apply the 
automated method to a wider sample of galaxies. 

We add these total hydrogen volume density measure- 
ments of M74 to the list of other spiral galaxies for which 
these measurements have been carried out, namely M33, 
M81, M83 and M101. With the automated pipeline de- 
scribed in this paper we will be able to do these measure- 
ments for a much larger sample of galaxies, using improved 
metallicity measurements as they become available. We also 
intend to compare these densities to other tracers of the in- 
terstellar medium as well as more global characterizations 
of galaxies in order to further validate the simple physical 
assumption of the presence of photodissociated atomic hy- 
drogen on galactic scales. If this assumption is indeed valid, 
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then these measurements provide an important independent 
check on CO molecular gas measurements, although addi- 
tional assumptions are needed to compare the two. 
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Figure 4. Total hydrogen volume densities across the disk of M74 near candidate PDRs as identified by clumpfind plotted on top of 
the GALEX far-UV image. The radius of the circles corresponds to the rounded logarithmic value of ntot- 




Figure 10. Left panel: Total hydrogen volume densities across the disk of M83 near candidate PDRs as identified by clumpfind plotted 
on top of the GALEX far-UV image. The radius of the circles corresponds to the rounded logarithmic value of n. Right panel: Locations 
of candidate PDRs in M83; Fig. 1 in H08b reproduced. 



Parameter Value 



Name NGC 628 (M74) NGC 5236 (M83) 

R.A. (J2000.0) 01:36:41.77 13:37:00.92 

Dec. (J2000.0) +15:47:00.5 -29:51:56.7 

Classification SA(s)c SAB(s)c 

Position angle 20° 43° 

Declination 7° 24° 

Distance 7.3 Mpc 4.5 Mpc 

R25 11.0 kpc 7.63 kpc 

E(B-V) 0.070 mag 0.066 mag 

Table 1. Basic properties of M74 and M83. Distance, position 
angle and declination of M74 as adopted in Warren c t al.| (|2010[) . 
The M83 parameters are identical to the ones adopted in[Hcincr 
et al.| (|2008b||. The M74 o ptical radius R25 is 5.2 arcmin as in 



Kamphuis &: Briggs|||1992||. Galac tic foreground extinctions E(B- 
V) are taken from |Schlegel et al.| ( |1998} . 
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Figure 5. Candidate PDR in M74. On the left hand side, con- 
tours mark the far-UV flux. On the right hand side, contours mark 
the HI column density. Both panels have the HI column density 
plotted in grayscale. The center of the region is marked with a 
star. Selected HI patches are marked with a white plus-sign. Fig- 
ures for all M74 candidate PDRs are available in the electronic 
version of this journal. 



Figure 1. Map of the 12 + log(0/H) values of the HII regions 
obtained from the SpIOMM observations. A color version of this 
map is available in the electronic version of this journal. 
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Figure 6. From top to bottom: M74 measured far-UV fluxes, HI 
column densities, separation between the central UV source and 
the candidate PDRs and incident fluxes Go- Fluxes are in units 
of ergs cm~ 2 s _1 A -1 . 



Figure 2. Sample interferogram of an HII region in M74. 
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Figure 3. Sample HII region spectrum showing the Ho and 
[Nil] lines as well as the SII doublet (not used here). The sky- 
subtracted spectrum is offset for clarity. 
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Figure 7. Dust-to-gas ratios at each candidate PDR. The 
black circles are the values taken from the SpIOMM mctallic- 
ity map, the values computed from a fixed metallicity slope are 
gray squares and the values taken from the Rosalcs-Ortega ct al. 
(2011) metallicity map are represented by gray crosses. 
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Figure 8. Total hydrogen volume densities for the three different 
metallicity scenarios. The symbols are the same as in Fig. [7] 
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Figure 9. Total hydrogen densities plotted against incident UV 
flux Go- The different selection effects are shown and described 
in the text. 



Figure 11. The measurements and results for M83. From top to 
bottom: dust-to-gas ratio, far-UV flux (in units of ergs cm -2 s — 1 
A" 1 ), HI column density, separation pni and incident flux Go- 
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Figure 12. The selection effects affecting the derived M83 total 
hydrogen volume densities. The individual effects are marked and 
explained in the text. 



© 2012 RAS, MNRAS 000,[T]j9] 



Hydrogen volume densities in nearby galaxies I 13 



4.0 



■ * 



■ • Bffim 

4 s Ob 5.53 S ■ ■ 

1° .1SW.!W9 safSg 



»« 9 ftC o V ^ o 

W™ □ ° B °«° 
•eg 1 □ a o °o 

a ■ ° 



0.1 



0.4 

R/R 2i 



Figure 13. Comparison of M83 total hydrogen volume densities 
from H08b and this work. The gray circles are the H08b densities. 
Open circles indicate densities obtained from a source contrast 
lower than 1. The results from this work (black circles) are only 
those with a source contrast above 1. 
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